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In the context of "two-metric" theories of gravity there is 
the possibility that cosmic defects will produce a spectrum 
of primordial adiabatic density perturbations. This will hap- 
pen when the speed characterising the defect-producing scalar 
field is much larger than the speed characterising gravity and 
all standard model particles. This model will exactly mimic 
the standard predictions of inflationary models, with the ex- 
ception of a small non-Gaussian signal which could be de- 
tected by future experiments. We briefly discuss defect evo- 
lution in these scenarios and analyze their cosmological con- 
sequences. 

PACS number(s): 98.80.Cq, 04.50.+h, 98.65.Dx, 98.70.Vc 



I. INTRODUCTION 

Cosmology is entering a crucial stage, where a growing 
body of high-precision data will allow us to determine a 
number of cosmological parameters, and to identify the 
mechanism that produced the "seeds" for the structures 
we observe today [Q. There are currently two classes of 
models that could be responsible for these — topological 
defect H and inflationary ||] models. The main differ- 
ence between them is related to causality. Initial condi- 
tions for the defect network are set up on a Cauchy sur- 
face that is part of the standard history of the universe. 
Hence, there will not be any correlations between quanti- 
ties defined at any two spacetime points whose backward 
light cones do not intersect on that surface. Inflation 
pushes this surface to much earlier times, and if the in- 
flationary epoch is long enough there will be essentially 
no causality constraints. This can also be seen by noting 
that inflation can be denned as an epoch when the comov- 
ing Hubble length decreases. It starts out very large, and 
perturbations can be generated causally. Then inflation 
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forces this length to decrease enough so that, even though 
it grows again after inflation ends, it's never as large (by 
today) as the pre-inflationary era value. Once primor- 
dial fluctuations are produced they can simply freeze in 
comoving coordinates and let the Hubble length shrink 
and then (for small enough scales) grow past them. 

As a step towards identifying the specific model that 
operated in the early universe, one would like to deter- 
mine which of the two mechanisms above was involved. 
The presence of super-horizon perturbations might seem 
a good enough test, but this is not the case: in defect 
models (as well as open or A-models) significant contri- 
butions are generated after the epoch of last scattering 
due to the integrated Sachs- Wolfe effect. The presence of 
the 'Doppler peaks' on small angular scales H is also not 
ideal: Turok || has shown that a causal scaling source 
can be constructed so as to mimic inflation and reproduce 
its contribution to the CMB anisotropies. This source 
is constructed "by hand" , and there is no attempt to 
provide a framework in which it could be realized. In 
any case, it shows that inflationary predictions are not 
as unique as one might think. We should also mention, 
however, that a nice argument due to Liddle Q (see also 
0) shows that the existence of adiabatic perturbations 
on scales much larger than the Hubble radius implies 
that either inflation occurred in the past, the pertur- 
bations were there as initial conditions, or causality (or 
Lorentz invariance) is violated. On the other hand, it is 
also possible to construct "designer inflation" models || 
that would have no secondary Doppler peaks, although 
these suffer from analogous caveats and they would still 
be identifiable by other means |j],|lo| . 

Finally, there are Gaussianity tests. There have been 
recent claims of a non-Gaussian component in the CMB 
flTf (but see also [0). Defects will generally produce 
non-Gaussian fluctuations on small enough scales [ p^[ , 
whereas the simplest inflationary models produce Gaus- 
sian ones. It's possible to build inflationary models that 
produce, eg. non-Gaussianity with a chi-squared distri- 
bution |Q, but if one found non-Gaussianity in the form 
of line discontinuities, then it is hard to see how cosmic 
strings could fail to be involved. 

This discussion shows that although defect and infla- 
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tionary models have of course a number of distinguishing 
characteristics, there is a greater overlap between them 
than most people would care to admit. It is also easy to 
obtain models where both defects and inflation generate 
density fluctuations (lj]Jl(| ■ The aim of this letter is to 
present a further example of this overlap. We discuss a 
model where the primordial fluctuations are generated by 
a defect network, but are nevertheless very similar to a 
standard inflationary model. The only difference between 
these models and the standard inflationary scenario will 
be a small non-Gaussian component. A detailed discus- 
sion will be presented in a forthcoming publication |0 . 



II. THE MODEL 

Our model follows the recent work on so-called 'vary- 
ing speed of light' theories jl^-Q, and more particularly 
the spirit of 'two-metric' theories |18 2^^,Q, having 
two natural speed parameters, say and c; the first is 
relevant for the dynamics of the scalar field which will 
produce topological defects, while the second is the or- 
dinary speed of light that is relevant for gravity and all 
standard model interactions. 

We assume that ^> c so that the correlation length 
of the network of topological defects will be much greater 
than the horizon size. We could, in analogy with p4] , p3[ , 
define our effective theory by means of an action, and 
postulate a relation between the two metrics. However, 
this is not needed for the basic point we're discussing, 
so we leave it for a future publication |l7[ ]. We concen- 
trate on the case of cosmic strings, whose dynamics and 
evolution are better known than those of other defects 
[^,^5-27 1 although much of what we will discuss will ap- 
ply to them as well. Note that could either be a con- 
stant (say [g^oo = ( c 0/ c2 )<7oo) or, as in [ p3| one could set 
up a model such that the two speeds are equal at very 
early and at recent times, and between these two epochs 
there is a period, limited by two phase transitions, where 
^S> c. As will become clear below, the basic mecha- 
nism will work in both cases, although the observational 
constraints on it will be different for each specific real- 
ization. 

The string network evolution is qualitatively analo- 
gous to the standard case [p|,p5|-[27] , and in particular a 
"scaling" solution will be reached after a relatively short 
transient period. The long-string characteristic length 
(or "correlation length") L will evolve as L — 'yc^t, 
with 7 = 0(1), while the string RMS velocity will obey 
Vcf, = /3c^, with j3 < 1. Note, however, that there 
are some differences relative to the standard scenario. 
The first one is obvious: if 3> c, the string network 
will be outside the horizon, measured in the usual way. 
Hence these defects will induce fluctuations when they 
are well outside the horizon, thus avoiding causality con- 
straints Note that compensation now acts outside the 



'c^-horizon'. We expect the effect of gravitational back- 
reaction to be much stronger than in the standard case 
[^8|,^3|. The general effect of the back-reaction is to re- 
duce the scaling density and velocity of the network rel- 
ative to the standard value (28). Thus we should expect 
fewer defects per "c^-horizon" , than in the standard case. 
However, despite this strong back-reaction, strings will 
still move relativistically. It can be shown J28| that al- 
though back-reaction can slow strings down by a mea- 
surable amount, only friction forces |2^,|27[] can force the 
network into a strong non-relativistic regime. Thus we 
expect v$ to be somewhat lower than c^, but still larger 
than c. Only in the case of monopoles, which are point- 
like, one would expect the defect velocities to drop below 
c due to graviton radiation [p3| . This does not happen 
for extended objects, since their tension naturally tends 
to make the dynamics take place with a characteristic 
speed Ccf, j2{|. This is actually crucial: if the network was 
completely frozen while it was outside the horizon (as in 
standard scenarios fl5|| ) then no significant perturbations 
would be generated. 

A third important aspect is that the the symmetry 
breaking scale, say E, which produces the defects can be 
significantly lower than the GUT scale, as density per- 
turbations can grow for a longer time than usual. The 
earlier the defects are formed, the lighter they could be. 
Proper normalization of the model will produce a further 
constraint on E. Finally, in the case where is a time- 
varying quantity which only departs from c for a limited 
period, the defects will become frozen and start to fall in- 
side the horizon after the second phase transition. Here 
we require that the defects are sufficiently outside the 
horizon and are relativistic when density fluctuations in 
the observable scales are generated. This will introduce 
additional constraints on model parameters, notably on 
the epochs at which the phase transitions take place. 



III. COSMOLOGICAL CONSEQUENCES 

In the synchronous gauge, the linear evolution equa- 
tions for radiation and cold dark matter perturbations, 
5 r and 6 m , in a flat universe with zero cosmological con- 
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where O a p is the energy-momentum tensor of the exter- 
nal source, 0+ = Ooo + On, a is the scale factor, "eg" 
denotes the epoch of radiation-matter equality, and a dot 
represents a derivative with respect to conformal time. 
We will consider the growth of super-horizon perturba- 
tions with cki] <C I. Then eqn. (Ill) becomes: 
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L + a -S m -U a -) 2 ( *t±*>*) S m = 4.G9+ , (3) 
a I \ a / \ a + a eq J 

and S r = 4(5 m /3. Its solution, with initial conditions 
S m = 0, S rn = can be written as 



^(x,7j) = 4ttG y */ I d 3 x'g(X;r,,r ) ')Q+(^,v'), (4) 



sin fcX , o ,, 
S( 5 77 , 77 )_ fc^ (5) 



Here X = |x — x'| and 'S' indicates that these are the 
'subsequent' fluctuations, according to the notation of 
p0{ , to be distinguished from 'initial' ones. 

We are interested in computing the inhomogeneities at 
late times in the matter era. When t/q ^> r] eq , the Green 
functions are dominated by the growing mode, oc ao/a eq , 
so the function we would like to solve for is IBOll 



T(k;r])= Urn —G(k, r/ , v) ■ 
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Consider the growth of super-horizon perturbations, for 
which the transfer function can be written |3(j| 
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Linear perturbations induced by defects are the sum of 
initial and subsequent perturbations: 

S m (k; Vo) = Vo) + s m( k ' Vo) 

fno _ 
= 4nG{l + z eq ) d V T c (k;v)e>+(k;v) , (8) 

Jr\i 

where is the epoch of defect formation. The trans- 
fer function for the subsequent perturbations, those gen- 
erated actively, was obtained in eqn. (Q) for super- 
horizon perturbations with ckr/Q <c 1. To include com- 
pensation for the initial perturbations, 5^, we make the 

substitution T c {k;rj) = (l + (fc c /fc) 2 ) T(k;rj), where 

k c oc (c^r])^ 1 is a long- wavelength cut-off at the com- 
pensation scale. This results from the fact that defect 
perturbations cannot propagate with a velocity greater 
than C0. For (c^o)" 1 <C k <C {c^rji)^ 1 the analytic ex- 
pression for the power spectrum of density perturbations 
induced by defects is 
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where !F(k, rj) is the structure function which can be ob- 
tained directly from the unequal time correlators f{0| -[32| . 
It can be shown |3l| that for a scaling network F{k, rj) — 
J-(krj) which, combined with the above relations gives 



P(k) oc / dr)S(krf)/r) oc k 
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where the function S, is just the structure function, J 7 , 
times the compensation cut-off function. Up until now 
we only considered the spectrum of primordial (ie, gener- 
ated at very early times) fluctuations induced by cosmic 
defects. In our model a Harrison-ZePdovich spectrum is 
predicted just as in the simplest inflationary models. The 
final processed spectrum will also be the same as for the 
simplest inflationary models. 

We investigate the Gaussianity of the string-induced 
fluctuations as in |l3j. The conclusions can easily be 
extended for other defect models. In the standard cos- 
mic string scenario the structure function T(k, rj) has a 
turn-over scale at the network correlation length, k^ = 
20 (c^rj)^ 1 p3| , |3l"| ] . At a particular time, perturbations 
induced on scales larger than the correlation length are 
generated by many string elements and are expected to 
have a nearly Gaussian. On the other hand, perturba- 
tions induced on smaller scales are very non-Gaussian 
because they can be either very large within the regions 
where a string has passed by or else very small outside 
these. This allows us to roughly divide the power spec- 
trum of cosmic-string-seeded density perturbations into 
a nearly Gaussian component generated when the string 
correlation length was smaller than the scale under con- 
sideration, and a strongly skewed non-Gaussian compo- 
nent generated when the string correlation length was 
larger (we call these the 'Gaussian' and 'non-Gaussian' 
contributions respectively). The ratio of this two com- 
ponents may be easily computed by splitting the struc- 
ture function in (g), in two parts: a Gaussian part 
!F g {k,rj) — T{k,rj) for k < k^ {T g = for k > k%) 
and a non-Gaussian part .F n g(fc,?7) = T(k,rj) for k > 
O^ng = for fc < fej). We can then integrate @ with 
this Gaussian/non-Gaussian split, to compute the rela- 
tive contributions to the total power spectrum. The final 
result will depend on the choice of compensation scale 
k c . If we take the maximum compensation scale allowed 
by causality J3J] (k c ~ 2 (c^rf)^ 1 ) the Gaussian contribu- 
tion to the total power spectrum will be less than 5%. In 
any case, the non-Gaussian contribution will always be 
smaller that the Gaussian one if, as expected, the com- 
pensation scale is larger or equal to the correlation length 
of the string network (k c < k(j). Departures from a Gaus- 
sian distribution are scale independent and analogous to 
those of standard defect models on large scales. 

By allowing for a characteristic velocity for the scalar 
field C0 much larger than the velocity of light (and grav- 
ity), we were able to construct a model with primordial, 
adiabatic (S r = 4<5 m /3), nearly Gaussian fluctuations 
whose primordial spectrum is of the Harrison-ZePdovich 
form. This is almost indistinguishable from the sim- 
plest inflationary models (as far as structure formation 
is concerned) except for the small non-Gaussian compo- 
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nent which could be detected with future CMB experi- 
ments. The Ci spectrum and the polarization curves of 
the CMBR predicted by this model should also be iden- 
tical to the ones predicted in the simplest inflationary 
models as the perturbations in the CMB are not gener- 
ated 'directly' by the defects. 



IV. DISCUSSION AND CONCLUSIONS 

We presented further evidence of the non-negligible 
overlap between topological defect and inflationary struc- 
ture formation models. The key ingredient is having the 
speed of the defect-producing scalar field much larger 
than the speed of gravity and standard model particles. 
This provides a 'violation of causality', as required by 
The only distinguishing characteristic of this model, by 
comparison with the simplest inflationary models, will be 
a small non- Gaussian signal. 

Admittedly our model could be considered "unnatu- 
ral" in the context of our present theoretical prejudices, 
and the same can certainly be said about other examples 
such as "mimic inflation" |j| and "designer inflation" || . 
Be that as it may, however, the fact that these examples 
can be constructed (and one wonders how many more are 
possible) highlights the fact that extracting robust pre- 
dictions from cosmological observations is a much more 
difficult and subtle task than many experimentalists (and 
theorists) believe. 



ACKNOWLEDGMENTS 

We thank Paul Shellard for useful discussions and 
comments. CM. is funded by FCT (Portugal) un- 
der 'Programa PRAXIS XXI', grant no. PRAXIS 
XXI/BPD/ 11 769/97. We thank CAUP for the facilities 
provided. 



[1] E. W. Kolb & M. S. Turner, The Early Universe, 
Addison- Wesley (1994). 

[2] A. Vilenkin & E. P. S. Shellard, Cosmic Strings and other 
Topological Defects, (Cambridge University Press: Cam- 
bridge, 1994). 

[3] A.D. Linde, Particle Physics and Inflationary Cosmology 

(Harwood, Chur, Switzerland, 1990). 
[4] A. Albrecht et al, Phys. Rev. Lett. 76, 1413 (1996); J. 

Magueijo et al, Phys. Rev. Lett. 76, 2617 (1996); W. Hu 

and M. White, Phys. Rev. Lett. 77, 1687 (1996). 
[5] N. Turok, Phys. Rev. D54, 3686 (1996); N. Turok, Phys. 

Rev. Lett. 77, 4138 (1996). 
[6] A. Liddle, Phys. Rev. 51, 5347 (1995). 



[-: 

[8 
[9 

[10 
[11 

[12 
[13 
[14 

[15 
[16 



[17 
[18 
[19 



[20 

[21 
[22 

[23 

[24 



[25 
[26 
[27 
[28 

[29 

[30 
[31 

[32 

[3:5 

[34 



Y. Hu, M. S. Turner and E. J. Weinberg, Phys. Rev. 
D49, 3830 (1994). 

D. Salopek, J. R. Bond and J. M. Bardeen, Phys. Rev. 
D40, 1753 (1989). 



A. Liddle, |astro-ph/9910110| (1999); J. Barrow and A. 
Liddle, Gen. Rel. Grav. 29 1501 (1997). 
A. Lewin and A. Albrecht, |astro-ph/990806l| (1999). 
P. Ferreira, J. Magueijo and K. M. Gorski, Ap. J. 503, 
LI (1998); J. Pando, D. Valls-Gabaud and L.-Z. Fang, 
Phys. Rev. Lett. 81, 4568 (1998). 

A. J. Banday, S. Zaroubi and K. M. Gorski, astro- 



ph/9908070 (1999) 



P. P. Avelino, E. P. S. Shellard, J. H. P. Wu and B. Allen, 
Ap. J. Lett. 507, L101 (1998). 

A. D. Linde, Phys. Lett. B158, 375 (1985); L. Kof- 
man, Phys. Lett. B173, 400 (1986); A. D. Linde and 
V. Mukhanov, Phys. Rev. D56, 535 (1997). 
P. P. Avelino, R. R. Caldwell and C. J. A. P. Martins, 
Phys. Rev. D59, 123509, (1999). 

C. Contaldi, M. Hindmarsh and J. Magueijo, Phys. Rev. 
Lett. 82, 2034 (1999); R. A. Battye and J. Weller, phys. 
Rev. D61, 043501 (2000). 

P. P. Avelino and C. J. A. P. Martins, in preparation. 

J. W. Moffat, Int. J. Mod. Phys. D2, 351 (1992). 

A. Albrecht and J. Magueijo, Phys. Rev. D59, 043516 

(1999); J. D. Barrow, Phys. Rev. D59, 043515 (1999); 

J. D. Barrow and J. Magueijo, Phys. Lett. B443, 104 

(1998). 

P. P. Avelino and C. J. A. P. Martins, Phys. Lett. B459, 
468 (1999). 



I. T. Drummond, |gr-qc/9908058| (1999). 

P. P. Avelino, C. J. A. P. Martins and G. Rocha, Phys. 

Lett. B483, 210 (2000). 

B. A. Bassett, S. Liberati, C. Molina-Paris and M. Visser, 



astro-ph/0001441 (2000) 



M. A. Clayton and J. W. Moffat, Phys. Lett. B460, 263 

(1999) ; M. A. Clayton and J. W. Moffat, Phys. Lett. 
B477, 269 (2 000); M. A. Clayton a nd J. W. Moffat 
gr-qc/OOO3O70| (2000); J. W. Moffat, |astro-ph/0003444 

(2000) . 

C.J. A. P. Martins and E.P.S. Shellard, Phys. Rev. D53, 
575 (1996). 

C.J. A. P. Martins and E.P.S. Shellard, Phys. Rev. D54, 
2535 (1996). 

C.J. A. P. Martins, Quantitative String Evolution, Ph.D. 
Thesis, University of Cambridge (1997). 
C. J. A. P. Martins, E. P. S. Shellard and B. Allen, 'The 
Velocity-dependent One-scale String Evolution Model', 
in preparation (2000). 

B. Carter, in Tlaxcala lecture notes, 2nd Mexican School 
on Gravitation and Mathematical Physics, A. Garcia et 
al. (eds.) (1996). 

S. Veeraraghavan and A. Stebbins, Ap. J. 365, 37 (1990) 
J.-H. P. Wu, P. P Avelino, E. P. S. Shellard and B. Allen, 
astro- ph/9812156| (1998) 



A. Albrecht and A. Stebbins, Phys. Rev. Lett. 68, 2121 
(1992). 

P. P. Avelino, E. P. S. Shellard, J.-H. P. Wu and B. Allen, 
Phys. Rev. Lett. 81, 2008 (1998). 

J. Robinson and B. D. Wandelt, Phys. Rev. D53, 618 
(1996). 



4 



